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We present strong numerical evidence for the band Kondo effect in the doped 2D Hubbard model by 
showing the systematic change ol the density of states, imaginary part of the self-energy, effective 
magnetic moment and quasiparticle residue upon decreasing the temperature. Quadratically vanish- 
ing (in frequency) scattering rates near the Fermi energy and a linearly vanishing (in temperature) 
effective magnetic moment at low temperatures strongly support the screening of magnetic moments 
from singly occupied electrons by doped holes in the 2D Hubbard model. 

PACS numbers: 71.10.Fd, 71.27.+a, 72.15.Qm 
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The Kondo effect [gjg] is realized by screening a lo- 
cal magnetic moment by conduction electrons below a 
characteristic temperature Tk- Below this temperature 
conduction electrons form a singlet with a localized elec- 
tron at the impurity site and as a result quench effectively 
magnetic degrees of freedom for a magnetic impurity. Re- 
cently a possibility of the so called band Kondo effect in 
the finite dimensional Hubbard model has received con- 
siderable attention, since its discovery in the infinite di- 
mensional Hubbard model [[| ||. In the latter limit, itin- 
erant electrons or doped holes play a role of conduction 
electrons and the other singly occupied electrons (due to 
a large Coulomb repulsion U) a role of local magnetic 
impurities. Consequently a sharp peak sensitive to the 
temperature appears at the Fermi energy in the density of 
states. In contrast to the usual Kondo or Anderson mod- 
els ||f7j , there is no clear distinction between conduction 
electrons in the Fermi sea and localized electrons at the 
impurity sites in this situation. The connection of this 
effect in the infinite dimensional Hubbard model with a 
Kondo-like screening is evident in that in this limit the 
model is reduced to an effective single impurity Anderson 
model where a Kondo or Abrikosov-Suhl resonance has 
been well established. Because of the absence of a formal 
connection in the finite dimensional model, however, it is 
quite uncertain whether the band Kondo effect survives 
or not in this case. 

Recently this problem was studied by Bulut et al. || by 
means of quantum Monte Carlo (QMC) calculations for 
the 2D Hubbard model. Although these authors showed 
a developing peak near the Fermi energy in the den- 
sity of states by decreasing the temperature, it is not 
clear whether this peak is associated with the onset of 
the Fermi liquid behavior or not, and more importantly 
whether it is due to the screening of magnetic moments 
(from singly occupied electrons) by doped holes. In this 
Letter, we will present strong numerical evidence for the 
band Kondo effect in the doped 2D Hubbard model by 
showing the systematic change of the density of states, 



imaginary part of the self-energy, effective magnetic mo- 
ment and quasiparticle residue upon decreasing the tem- 
perature. 

There are fundamental differences between the two and 
infinite dimensional Hubbard models. In contrast to the 
latter model where an Abrikosov-Suhl resonance appears 
for any infinitcsimally small doping of a half-filled band, 
strong 2D spin fluctuations near half-filling Hfio) easily 
destroy the singlet formation in the 2D model. Hence, 
forming a quasiparticle state at the Fermi energy, if it is 
possible, is expected to happen only below a certain elec- 
tron concentration in the 2D Hubbard model. In order 
to take into account strong 2D critical spin fluctuations 
properly, we impose the following three exact sumrules 
to the spin, charge, and particle-particle susceptibilities 
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T and N are the absolute temperature and number of 
lattice sites, q is a compact notation for (q,iv n ) where 
iv n are either Fermionic or Bosonic Matsubara frequen- 
cies. The dynamical spin, charge and particle-particle 
susceptibilities are calculated by 
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Xph(q) an d Xp P (q) are irreducible particle-hole and 
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particle-particle susceptibilities, respectively, which are 
computed from 

&(9) = -J£g°(*-?)g°(*) 
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where G (k) is the noninteracting Green's function. 
U sp , U c h, and U pp in Eq. ^ are renormalized inter- 
action constants for each channel which are calcu- 
lated self-consistentlyby making an ansatz U sp = 
(7(n T n i )/((n T )(n i )) @ in Eq. | By defining U sp , U ch , 
and f/pp this way, the Mermin- Wagner theorem |ll[] as 
well as correct atomic limit for large u> are satisfied simul- 
taneously (|] . In order to find the chemical potential for 
interacting electrons, first we calculate Eqs. (0)-(||) and 
the self-energy (Eqs. (Q) in Ref. ||) with the noninteract- 
ing Green's function whose noninteracting chemical po- 
tential gives a desired electron concentration. Then, the 
chemical potential for interacting electrons is determined 
in such a way that the calculated electron concentration 
with the interacting Green's function is the same as the 
desired value. Throughout the calculations the unit of 
energy is t and all energies are measured from the chemi- 
cal potential \i. We used a 128 x 128 lattice in momentum 
space and performed the calculations by means of well- 
established fast Fourier transforms (FFT). It should be 
also noted that we used a real frequency formulation in 
Eqs. (|lj)-(|j) to avoid any possible uncertainties associated 
with numerical analytical continuation. 

We start in Fig. [I] by studying the density of states 
at various electron concentrations n = 1.0, 0.91, 0.85 
and 0.80 (Fig. 13(a)) and at n — 0.80 for various tem- 
peratures (Fig. Wb)) for [7 = 8. Within our numerical 
calculations, we do not have the self-consistent solution 
for U sp at lower temperatures than 0.05 near half-filling. 
Note that the spin-spin correlation length exponentially 
increases like £ ~ exp(constant /T) at low temperatures 
near half-filling. At n = 1.0 (dot-dashed curve), strong 
2D critical spin fluctuations due to the Mermin- Wagner 
theorem deplete the spectral weight near the Fermi en- 
ergy, leading to an antiferromagnetic pseudogap. Since 
the imaginary part of the self-energy is proportional to 
£ at u> — for low temperatures near half-filling 
scattering rates increase exponentially with decreasing 
temperature, eventually leading to an antiferromagnetic 
insulator at zero temperature. At n = 0.91 (dashed 
curve) , strong spin fluctuations are still dominating near 
the Fermi energy as the strong persistence of the antifer- 
romagnetic pseudogap indicates. Decreasing the temper- 
ature still gives rise to exponentially increasing scattering 
rates. This is in clear contrast to the infinite dimensional 
model where a sharp quasiparticle peak associated with 
an Abrikosov-Suhl resonance appears by infinitesimally 
small doping of a half-filled band. At n = 0.85 (dotted 
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FIG. 1. Density of states for U = 8 (a) at n = 1.0, 0.91, 
0.85 and 0.80 denoted as the dot-dashed, dashed, dotted, and 
solid curves, respectively and (b) at T = 0.8, 0.4, 0.2, 0.1 and 
0.05 denoted as the dot-dashed, long-dashed, dashed, dotted, 
and solid curves, respectively at n = 0.80. In (a) for n = 1.0 
and 0.91, T is 0.05, while for n = 0.85 and 0.80, T is 0.005. 

curve), most of the antiferromagnetic pseudogap is closed 
and a small peak develops at the Fermi energy. A care- 
ful investigation of the imaginary part of the self-energy 
(not shown in this Letter), however, still shows a small 
remnant of the 2D critical fluctuations near the Fermi 
energy. At n — 0.80 (solid curve), the antiferromagnetic 
pseudogap is completely closed by doping and a narrow, 
sharp spectral weight appears at the Fermi energy. 

In Fig. 0(b), the systematic change of the density of 
states for U — 8 and n = 0.80 is presented for T = 0.8 to 
0.05. As the temperature is decreased, the general shape 
of the density of states remains unchanged except near 
the Fermi energy. A small spectral weight in the interme- 
diate frequency regime is transferred near the Fermi en- 
ergy and a featureless background at high temperatures 
is transformed into a sharp peak at the Fermi energy 
with decreasing temperature. In order to find that this 
feature is associated with the onset of the Fermi liquid, 
we present the imaginary part of the self-energy in Fig. ^ 
for the same parameters as in Fig. |l|(b) . 

First ImYi{k,uj) is shown in the wide range of fre- 
quency axis in Fig. ^(a). With decreasing temperature 
scattering rates are progressively decreased particularly 
near the Fermi energy. The low frequency behavior of 
the scattering rates is shown in Fig. ||(b). Below T — 
0.1 — 0.2, the imaginary part of the self-energy vanishes 
quadratically near the Fermi energy with a small con- 
stant shift due to a finite temperature effect. The same 
quantity for T = 0.005 is also plotted as the dot-dashed 
curve where the quadratically vanishing scattering rates 
are more visible. The log-log plot of the scattering rates 
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FIG. 2. Imaginary part of the self-energy for U — 8 and 
n = 0.80 at the Fermi surface (a) for T = 0.8, 0.4, 0.2, 0.1 and 
0.05 denoted as the dot-dashed, long-dashed, dashed, dotted, 
and solid curves, respectively in the wide range of frequency 
axis and (b) for T = 0.4, 0.2, 0.1 and 0.05 denoted as the 
long-dashed, dashed, dotted, and solid curves, respectively in 
the low frequency regime. The dot-dashed curve in (b) is for 
T = 0.005. 

vs. frequency for T = 0.005 shows JmS(fcj?, w) ~ u 1 
in an interval of ±[0.06 — 0.61]. 

In order to establish a more firm ground that this is in- 
deed due to the screening of magnetic moments by doped 
holes, we show the effective magnetic moment defined as 
Txsp(0, 0) upon decreasing the temperature (filled circles 
in Fig. ||(a)). Below T = 0.2, the effective magnetic mo- 
ment vanishes linearly in temperature, a clear indication 
of the Kondo screening. Above T = 0.2, the effective 
magnetic moment deviates significantly from a straight 
line and appears gradually saturating at high tempera- 
tures. The quasiparticle residue is also presented as the 
open circles in Fig. ^(a). With decreasing temperature 
it saturates at 0.1 around T = 0.05. Hence, the effective 
mass of the quasiparticle becomes ten times heavier than 
the bare electron mass below T = 0.05. 

We also performed the same calculations for U = 4. 
Due to a smaller coupling strength, the complete de- 
struction of the 2D critical spin fluctuations happens at 
n = 0.87. In Fig. |](a) the density of states for U = 4 
and n = 0.87 is shown at T = 0.8 to 0.05. The general 
shape of the density of states for U = 4 is quite different 
from that for U = 8 but instead closer to that of nonin- 
teracting electrons. With decreasing temperature some 
spectral weight in the intermediate frequency regime is 
transferred near the Fermi energy to build a sharp peak 
which can be interpreted as a Kondo resonance. We 
present the scattering rates in the low frequency regime 
in Fig. |(b). As in Fig. |(b), below T = 0.1 - 0.2, the 
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FIG. 3. Effective magnetic moment (filled circles) and 
quasiparticle residue (open circles) (a) for U = 8, n = 0.80 
and (b) for U = 4, n = 0.87. The solid lines are a linear 
interpolation of the effective magnetic moment at low tem- 
peratures. 

imaginary part of the self-energy vanishes quadratically 
near the Fermi energy with a small constant shift. The 
same quantity for T — 0.005 is also presented as the dot- 
dashed curve. The log-log plot of the scattering rates vs. 
frequency for T = 0.005 shows JmS(fei?, u>) ~ to 1 - 95 in an 
interval of ±[0.05 - 0.37]. In Fig. |(b), we show the ef- 
fective magnetic moment (filled circles) and quasiparticle 
residue (open circles) for U — 4. The effective magnetic 
moment for {7 = 4 also vanishes linearly in temperature 
below T = 0.2. The quasiparticle residue saturates at 
0.29 below T = 0.025. 

Finally we show in Fig. || the density of states and 
imaginary part of the self-energy for U = 1 and n = 0.95 
where critical spin fluctuations are completely destroyed. 
For this weak coupling strength, magnetic moments from 
singlely occupied electrons are not well developed and 
thus the screening of the magnetic moments by doped 
holes is questionable. The density of states from T = 0.8 
to 0.05 in Fig. |](a) is very close to that for noninteract- 
ing electrons with a small constant shift in the frequency 
axis. Within our numerical resolution it is difficult to find 
any developing peak near the Fermi energy by decreasing 
the temperature. The imaginary part of the self-energy 
in Fig. ||(b), however, clearly exhibits the drastic change 
of the scattering rates to the quadratic behavior near the 
Fermi energy with decreasing temperature. The log-log 
plot of the scattering rates vs. frequency for T — 0.005 
shows ~ uj 1 - 97 in an interval of ±[0.04—0.18]. 

The quasiparticle residue is found to saturate at 0.92 and 
the effective magnetic moment also vanishes linearly in 
temperature with decreasing temperature. All these find- 
ings for U = 1 are also consistent with the existence of 
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FIG. 4. (a) Density of states for T = 0.8, 0.4, 0.2, 0.1 and 
0.05 denoted as the dot-dashed, long-dashed, dashed, dotted, 
and solid curves, respectively and (b) imaginary part of the 
self-energy at the Fermi surface for T = 0.4, 0.2, 0.1 and 
0.05 denoted as the long-dashed, dashed, dotted, and solid 
curves, respectively in the low frequency regime for U — 4 
and n = 0.87. The dot-dashed curve in (b) is for T = 0.005. 



FIG. 5. (a) Density of states for T = 0.8, 0.4, 0.2, 0.1 and 
0.05 denoted as the dot-dashed, long-dashed, dashed, dotted, 
and solid curves, respectively and (b) imaginary part of the 
self-energy at the Fermi surface for T = 0.4, 0.2, 0.1 and 
0.05 denoted as the long-dashed, dashed, dotted, and solid 
curves, respectively in the low frequency regime for U = 1 
and n = 0.95. The dot-dashed curve in (b) is for T = 0.005. 



the Kondo screening in the doped 2D Hubbard model 
found earlier for larger interaction strengths. 

In summary, we have presented strong numerical evi- 
dence for the band Kondo effect in the doped 2D Hub- 
bard model by showinging the systematic change of the 
density of states, imaginary part of the self-energy, ef- 
fective magnetic moment and quasiparticle residue upon 
decreasing the temperature. Quadratically vanishing (in 
frequency) scattering rates near the Fermi energy and 
a linearly vanishing (in temperature) effective magnetic 
moment at low temperatures strongly support the screen- 
ing of the magnetic moments from singly occupied elec- 
trons by doped holes in the doped 2D Hubbard model. 
Our numerical calculations indicate that the band Kondo 
effect in the doped 2D Hubbard model persists in the 
strong to weak coupling regimes. 
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